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I.  INTRODUCTION 


Recently,  a considerable  interest  has  developed  in  the  temperature 
fluctuations  occurring  within  heated  boundary  layers.  The  specific  important 
applications  include  the  effects  of  such  temperature  fluctuations,  notably 
through  attendant  density  fluctuations,  upon  laser  beam  propagation  through 
boundary  layers.  Primarily,  such  problems  arise  for  beams  being  sent 
through  boundary  layers  occurring  about  jet  planes,  flying  at  considerable 
altitude,  and  at  speeds  near  (the  ambient)  Mach  1.  The  density  changes  that 
affect  laser  propagation  can  be  caused,  of  course,  either  by  the  temperature 
fluctuations,  considered  here,  or  by  pressure  fluctuations;  the  latter  effects 
have  been  treated  elsewhere^  by  techniques  similar  to  those  used  in  this 
paper.  We  assume  in  this  work  that  the  temperature  effects  dominate,  and, 
we  neglect  the  pressure  fluctuation  effects. 

Usually  work  on  heated  boundary  layers  has  been  directed  toward 
finding  rms  temperature  fluctuations  (among  other  quantities)  at  the  most. 
Here  for  the  laser  propagation  problem,  and  for  other  problems  as  well, 
the  temperature  correlation  function  is  needed.  For  simple  rms  values  of 
the  temperature  fluctuation  it  would  be  possible  to  use  some  variation  of 
mixing-length  methods.  But  of  course  for  the  correlation  function  much 
more  information  is  needed.  Our  approach  here  (involving  beginning  with 
the  usual  equations  of  motion  for  the  fluid)  uses  an  approximate  Ijagr2uigian 
integration  procedure  in  order  to  obtain  the  required  results. 
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The  basic  procedure  to  be  followed  consists  of  beginning  with  the 
energy  equation  for  fluid  flow.  In  that  equation  we  shall  genera  Ty  neglect 
the  effects  of  thermal  conduction  within  the  heated  turbulent  boundary  layer, 
with  the  exception  of  the  region  very  near  the  wall  bounding  the  layer;  in 
this  (laminar)  sublayer  and  regions  contiguous  to  it,  there  is  heat  generated 
by  the  viscous  action,  as  a result  of  the  very  high  velocity- shear  occurring 
very  near  the  wall.  This  heat  generation,  given  quantitatively  by  the  dissipa- 
tion function,  is  included  in  the  problem  as  part  of  the  heat  flux  occurring 
within  the  boundary  layer.  In  addition,  there  is  a rapid  deceleration  of  the 
fluid  near  the  wall  through  viscous  action;  since  this  process  is  riot  adiabatic, 
there  is  some  irreversible  heat  generation,  also  to  be  included  in  the  heat 
flux.  And,  finally,  in  problems  of  this  type  the  wall  temperature  is  in  general 
different  than  the  ambient  temperature.  In  the  application  cited  above,  it  is 
considerably  higher  than  the  external  atmospheric  temperatures.  In  our 
treatment,  this  heated  wall  is  also  taken  into  account.  This  effect  can,  again, 
be  included  as  part  of  the  heat  flux  within  the  boundary  layer. 

The  boundary  layer  near  a flat  plate  or  wall  is  initiated  at  the  leading 
edge  and  grows  slowly  as  one  progresses  downstream  (in  the  x-direction). 

In  our  treatment,  we  assume  (as  is  usual)  that  the  layer  growth  rate  is  slow. 
The  process  does  not  change  appreciably,  in  regions  of  interest,  in  the 
downstream  direction.  The  derivatives  of  important  quantities  within  the 
flow  field  are  much  greater  in  the  normal  direction  (the  z-direction)  than  in 
the  x-direction.  For  such  a flow,  the  velocity  fluctuations  a.e  approximately 
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statistically  homogeneous  in  planes  parallel  to  the  wall.  Further,  even  for 
flows  approaching  Mach  1,  it  is  well  known  that  the  velocity  fluctuations, 
which  are  particularly  important  in  the  treatment  here,  are  at  relatively 
small  Mach  number:  the  rms  velocity  is  considerably  below  the  mean  flow 
velocity,  typically  no  more  than  10%  of  that  velocity.  Thus,  for  the  fluctua- 
tions, one  can  treat  the  problem  as  incompressible.  We  carry  that  assump- 
tion here  but,  of  course,  allow  temperature  variations  to  produce  attendant 
density  variations.  In  effect,  we  are  adopting  the  familiar  assumption  that 
the  dynamics  of  the  process  are  but  slightly  influenced  by  compressible 
effects.  It  will  be  seen  in  the  development  that  the  temperature  fluctuations 
under  our  assumptions  can  be  determined  from  the  velocity  fluctuations  by 
an  integration  process.  We  do  not  attempt  here  to  solve  for  velocity 
fluctuations:  we  suppose  that  they  are  given. 

In  this  treatment,  we  assume  that  the  total  convected  heat  flux  is  a 
given  quantity  and  that  quantity  will  be  determined  from  familiar  fluid 
mechanics  theory. 

The  needed  velocity  correlations  are  taken  from  statistically 
homogeneous  and  isotropic  flows  (notably  wind  tunnel  experiments),  but  the 
correlations  are  weighted  with  the  rms  velocity  fluctuations  at  the  measuring 
points  for  the  velocity  field.  The  scale  for  the  correlation  functions  is 
proportional  to  the  distance  from  the  wall  so  long 'as  one  lies  within  the 
boundary  layer  and  is  6,  the  displacement  thickness  of  the  boundary  layer, 
outside  that  layer. 
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We  proceed  to  calculate  the  temperature  space- correlation  (and  its 
variance).  We  must  be  given  the  turbulence  intensity  level  in  this  treatment. 
A separate  calculation  (using  standard  methods  from  heat  transfer  theory) 
is  made  to  obtain  the  convective  heat  flux,  just  above  the  layer  6^. 


n.  THEORY 


A sketch  of  the  boundary  layer  problem  is  shown  in  Fig.  1,  where  the 
average  velocity  u depends  only  on  z,  following  our  parallel  flow  assumption 
for  the  turbulent  boundary  layer.  For  this  problem,  the  free  stream  velocity 
is  u and  6 is  the  displacement  thickness  of  the  boundary  layer  defined  by 


t 

L 


(1) 


The  average  temperature,  also  a function  of  z alone,  is  T(z).  The  tempera- 


ture of  the  wall  is  T . The  laminar  sublayer  is  of  thickness  6 . In  Fig.  1, 

'n  8 


the  heating  effect  of  the  dissipation  and  the  imperfect  temperature  recovery, 
confined  to  regions  near  the  sublayer,  is  shown  to  cause  an  increase  in  the 
average  temperature.  The  ambient  value  of  the  temperature,  T^,  is  shown 
for  z large.  The  temperature  T is  assumed  here  to  exceed  the  ambient 


value  though,  of  course,  if  it  were  less  than  that  value  the  treatment  would 


be  the  same.  The  energy  equation  for  compressible  fluid  flow  is  given  by 


(2) 


where  the  dissipation  function  9 is  given  by 


* = ^ [("x)  + + ('^z)  ] + (^x  + '^y)  + K + 


(3) 
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♦ ♦ ♦ 

u = (u  , V , w ) 


and  where  subscripts  x,  y,  and  z indicate  partial  derivatives  here  and  later. 
We  also  use  the  corresponding  component  notation  for  other  fields  later. 

For  incompressible  flow  Eq.  (2)  reduces  to  (for  approximately  constant 
thermal  conductivity  k) 


DT  2 ♦ 

pCp  ^ + p® 


where  ® is  given  by  Eq.  (3)  without  the  last  term.  In  these  equations,  p, 
assumed  here  to  be  approximately  constant,  is  the  density,  Cp  is  the  specific 
heat  at  constant  pressure,  per  unit  mass,  and  p is  the  viscosity  of  the  fluid. 
Further,  D/Dt  is  the  material  derivative  given  by 


DT*/Dt  = eT*/et  + u*.  VT* 


or,  for  incompressible  flow,  if  convenient,  we  can  write 


DT*/Dt  = 8T*/at+V  • (u*T*) 


The  change  in  density  is  assumed  small  in  our  work.  In  Eq.  (4),  strictly 
speaking,  for  incompressible  flow  there  is  but  one  specific  heat.  However, 
we  assume  here  that  there  is  a response  of  the  medium  to  the  temperature 
changes  set  in  near  the  wall.  The  appropriate  specific  heat  is  that  at 
constant  pressure,  for  such  a slight  variation  from  incompressible  flow. 
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The  dissipation  function  is,  as  suggested  earlier,  small  throughout 
most  of  the  boundary  layer  and  may  be  neglected  there.  In  the  region  6^ 
near  the  wall,  because  of  the  very  high  shear- velocity,  there  may  be  an 
appreciable  generation  of  heat.  Accordingly,  we  treat  the  region  just  outside 
the  viscous  sublayer  and  take  into  consideration  the  heat  generated  within 
that  layer  when  calculating  the  total  heat  transfer  through  the  boundary  layer. 
Furthermore,  we  neglect  (for  either  air  or  water)  the  very  slight  heat  con- 
duction term,  which  is  small  compared  with  the  convective  effects,  outside 
the  sublayer  In  dimensionless  form,  we  know  that  the  heat  conduction 

is  of  order  the  viscous  dissipation.  For  large  Reynolds  number  flows  being 
considered  here,  these  two  effects  are  in  general  small  except  for  special 
regimes:  specifically,  the  heat  conduction  is  of  order  the  convective  effect, 
near  the  laminar  sublayer  6^.  The  result  of  these  considerations  is  that  our 
governing  equation  becomes 

DT*/Dt  =0  (6) 

Thus,  there  is  no  change  in  temperature  (no  heat  exchange)  following  the  flow. 
Under  the  stated  assumptions,  Eq,  (6)  may  be  expected  to  apply  everywhere 
outside  fig  — through  the  boundary  layer  and  out  to  regions  beyond  the  layer. 
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m.  DERIVATION  OF  AVERAGES  AND  CORRELATION  FUNCTIONS 


We  begin  by  writing  the  field  quantities  in  terms  of  their  averages  and 
fluctuations  about  those  averages. 


T*  = T + T;  u*  = u + u 


with,  by  definition. 


<T>  = <u>  = 0 


By  incompressibility  we  see  that 


y.u  =y,u=y.u=0 


The  overbar  indicates  the  (time)  average;  we  assume  that  the  process  is 
statistically  stationary.  Substituting  Eq.  (7)  into  Eq.  (6)  and  averaging  the 
resulting  equation,  we  find  that 


flT/8t  + V • (uT  + uT)  = 0 


Note  that  the  heat  flux  q convected  in  the  normal  direction,  away  from 


the  wall,  is  given  by 


q/Cpp  = ^ 


One  can  see  this  relation  by  considering  the  followirs:  If  a parcel  of  fluid 
with  slightly  higher  temperature  is  tremsported  upward,  that  is  with  w 
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positive,  there  will  be  an  increase  in  heat  in  the  upper  region,  and 

i 

conversely.  Thus,  the  net  heat  flux  flowing  upward  is  proportional  to  wt. 


We  again  substitute  Eq.  (7)  into  Eq.  (6)  and  subtract  Eq.  (9)  to  find 
the  fluctuation  relation,  using  Eq.  (8). 


T + (u  +u).VT  = - V.  (uT  - uT) 


(li) 


Equation  (11)  states  that  the  material  derivative  of  the  temperature  fluctua- 
tion (the  time  rate  of  change  of  that  fluctuation  following  the  flow)  is  equal 
to  the  two  terms  on  the  right  side.  Of  those  terms  the  first  is  considerably 
larger  than  the  second,  involving  as  it  does  the  average  temperature  rather 
than  the  fluctuation;  we  neglect  the  second  term  here.  Before  attempting  to 
integrate  this  equation,  some  discussion  is  necessary.  Begin  with  a parcel 
of  fluid  at  constant  temperature,  that  is,  with  no  fluctuation;  such  constant 
temperature  boundaries  occur  at  the  wall  and  in  the  free  stream.  This 
parcel  of  fluid  then  diffuses  by  what  amounts  to  a random  walk  into  the 
boundary  layer,  arriving  at  the  point  of  interest.  This  diffusion  path  is 
nearly  horizontal  but  of  course  not  exactly  so.  The  vertical  distance 
through  which  it  diffuses  in  time  Tj^  is  of  order  Az, 


(Az)^  =<(/*"  > S (w')^  Tj5  t 


(12) 


where  t^  is  the  characteristic  time  for  the  turbulence  process  and  the  prime 
as  usual  indicates  the  rms  value  of  the  vertical  velocity  fluctuation. 


i 


-i4- 


The  characteristic  time  for  the  turbulent  process  is  known  to  be  approximately 

to  = 

The  distance  Az,  which  is  required  here  for  the  diffusion  process,  is  quite 
complicated.  Near  the  wall  the  distance  is  essentially  z,  the  distance  from 
the  wall,  and  the  fluid  which  arrives  consists  mainly  of  hot  elements  of 
fluid.  On  the  other  hand,  near  the  edge  of  the  boundary  layer,  that  is,  z 
near  8,  the  diffusion  distance  is  something  like  8-z  and  the  fluid  that  arrives 
comes  mainly  from  the  cold,  free  stream.  Finally,  no  fluid  with  temperature 
fluctuation  arrives  at  points  outside  the  boundary  layer,  so  tiiat  there  the 
diffusion  time  may  be  taken  as  zero.  We  use  for  the  distsmce  over  which  the 
diffusion  occurs,  the  convenient,  approximate  expression 

-2  -2  -2 
Az  = z + (8  - z)  , z < 8 

= 0 , z>  8.  (13) 

Combining  the  simple  expressions,  we  obtain  for  the  time  over  which  the 
diffusion  occurs,  the  result 

Tjj  = j[z"^  + (8  - z)"^]  8w'j"*,  z < 8 (14) 

and  Tj^  = 0 for  z > 8 . Using  the  incompressibility  condition  on  u and  the 
parallel  flow  assumption,  we  can  write  for  Eq.  (11),  neglecting  the  term 
- V • uT, 


DT/Dt  = -wT 

z 


(15) 


where  the  subscript  indicates  differentiation;  now,  we  formally  integrate 
that  equation  to  find 


T(r) 


- u(*)t,  y. 


z;  t - 


(16) 


We  suppose  that  the  parcel  of  fluid  began  in  a region  where  there  was  no 
temperature  fluctuation.  In  integrating  along  the  path  of  the  fluid  as  shown 
in  Eq.  (16)  we  have  dropped  the  velocity  fluctuation  in  the  argument  of  2 
since  it  in  general  is  much  smaller  than  the  mean  flow  used  there.  Dropping 
the  velocity  fluctuation  has  the  effect  of  integrating  the  function  w along  a 
line  parallel  to  the  x-axis  rather  than  following  the  actual  path  taken  by  the 
fluid  particle.  This,  of  course,  introduces  some  error  in  the  result  (see 
Appendix  A).  The  various  accumulated  errors  in  the  treatment  will  be 
handled  by  introducing  a multiplicative  constant,  which  will  be  adjusted  in 
order  to  give  the  proper  heat  flux  near  the  laminar  sublayer. 

We  proceed  to  construct  the  temperature  correlation  function.  To  do 
this,  we  write  Eq.  (16)  at  the  point  r'  rather  than  r as  here,  then  do  the 
same  at  the  point  r' '.  Multiplying  the  resulting  equations  together  and 
averaging,  we  obtain 


T(r',  r")  ■ <T(r')  T(r")> 


= T ,T 


T 


(17) 


X R33(x  - u(z')  (t  + t')  + u(z')t',  y;  z',  z",  t] 
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In  Eq.  (17)  we  used  the  parallel  flow  assumption,  that  the  velocity  correlation 
is  statistically  homogeneous  in  planes  parallel  to  the  wall.  As  a consequence, 
the  correlation  is  a function  of  the  difference  arguments  in  the  x~  and  y” 
directions.  That  correlation  is  defined  by  the  relation 


where  we  employ  the  notation 


r = r"  - r' 


We  also  changed  one  variable  of  time  integration  from  t^'  to  t = t''  - t'. 

For  Eq.  (17)  we  know  that  the  scale  of  the  correlation  function  in  space  is  of 
order  8,  and  in  time  is  of  order  t^,  given  by  Eq,  (12).  It  is  seen,  tiien, 
that  the  integral  is  controlled  by  the  integration  of  tiie  time  dependence  in 
the  x-argument.  The  time  dependence  in  tiie  time  argument  itself  has  a 
much  larger  scale.  Using  that  fact  we  can  integrate  Eq.  (17)  and  obtain  the 
result 


T(r',  r")  = T , T „ (I,  y;  z',  z")  dg  (20) 

- - * * y.«  33 

where  may  be  taken  to  be  the  simultaneous  (z-directed)  velocity  correla- 
tion, In  Eq,  (20),  there  is  no  x-dependence  for  ttie  temperature  correlation 
until  the  separation  between  the  measuring  points  in  the  x-direction  exceeds 
(the  erroneous  x-scale)  uTjj~  uM/w^  Thus,  tfie  scale  in  the  x-direction 
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would  be  considerably  greater  than  the  scale  in  the  other  directions.  This, 

as  is  pointed  out  in  Appendix  A,  is  a result  of  the  approximation  made  in 

Eq.  (16),  Further  as  pointed  out  there,  a multiplicative  constant  is  introduced 

2 

by  the  approximation,  inserted  in  Eq.  (20)  and  called  K . 

In  order  to  proceed  with  the  calculation,  we  need  the  correlation  R, 
we  make  the  plausible  assumption  that  it  is  approximately  given  by 


‘■33’ 


= w'(z")  yff'iz")  R33(r) 


(21) 


The  R^j  is  the  normalized,  second-order  velocity  correlation  for  incompres- 
sible, homogeneous  and  isotropic  turbulence.  These  correlations  have  been 
much  measured  experimentally.  It  is  seen  that  the  proposed  velocity  correla- 
tion for  this  problem  is  a weighted,  statistically  homogeneous  correlation, 
the  weighting  being  taken  at  the  measuring  points  for  the  correlation  function. 

There  is  a question  as  to  whether  an  incompressible  flow  field  may  be 

represented  using  such  weighting  functions.  It  is  shown  in  Appendix  B that 

this  is  possible.  By  our  parallel  flow  assumption  the  rms  velocities,  being 

average  quantities,  depend  only  on  z.  From  the  theory  of  homogeneous 

3 

and  isotropic  turbulence  we  know  that  the  normalized  correlation  can  be 
written 


.*  _ 1 ®^(m)  r r + 

ij  - -7¥-81-  Vj  + 


R..  = 


f + 1/2  r 


8f 

77 


hij 


(22) 
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^33  = ^(tI)  + (u)  f - ■^) 


(23) 


with 


f*r\m 

* \ W TJtJYX) 


where  M is  the  integral  scale  of  the  turbulence,  and  f is  the  longitudinal 
correlation.  Referring  again  to  Ref.  3,  the  transverse  correlation  is 
related,  to  the  longitudinal 


(24) 


Consider  M,  the  scale  of  the  turbulence.  Near  the  wall,  the  scale  is 
given  by  the  distance  from  the  wall.  For  the  two-point  correlations  con- 
sidered here  the  scale  must  be  given  by  the  lesser  distance.  As  we  proceed 
through  the  boundary  layer,  finally  the  scale  approaches  the  thickness  of  the 
boundary  layer  and  outside  that  layer  we  can  expect  that  any  disturbances 
which  remain  will  retain  the  scale  of  the  boundary  layer  thickness.  We 
accordingly  choose  for  the  function  M,  the  following  form 


M"^(*',  z")  = (*')"^  + (*")"^  + 6"^ 


(25) 
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Jit 


a 8^  (28) 

ThenT^/  6^  may  be  approximated  by  - T^.  We  obtain  for  the  mean  ^ 

square  temperature  fluctuation,  the  form 


h 


or,  using  the  nominal  value  =0.05  (Ref.  7) 


T'(ig)  = (T^  - TJ  K{0.05)^/^ 


Here,  a is  the  integral 


05 

a=f  f<§)d5=0.5 


the  numerical  value  for  this  integral  is  obtained  below  from  a large  Reynolds 
number  calculation  for  the  longitudinal  correlation  function  f.  As  suggested, 
we  compare  this  result  with  that  obtained  from  heat  transfer  theory  given  in 
Eq.  (Cl 3).  The  result  for  K is 


K = 0.  16  (M  x)‘ 


The  rms  temperature  fluctuation  throughout  the  layer  is. 


T'(z')  = 0.12  T ,(z 


from  Eq.  (27),  where  is  Im  and  x is  measured  from  the  beginning  of  the 
boundary  layer. 

For  simplicity,  we  assume  (without  great  error)  that  the  mean 
temperature,  the  turbulent  velocity  fluctuation,  and  the  mean  velocity  do  not 


r 


vary  greatly  in  going  from  z'  to  z",  then  the  normalized  temperature 
correlation  function,  found  by  dividing  Eq.  (27)  by  the  square  of  Eq.  (33), 
is 


[T 


^ I 


f(p)  + |,  f '(p) 


= t C^)  + + C^) 

where 


(34) 


I 
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IV.  DISCUSSION  OF  RESULTS  AND' CONCLUSIONS 


The  standard  deviation  of  the  temperature  (the  rms  fluctuation)  given  in 
Eq.  (33)  has  a constant  chosen  to  match  the  expected  standard  deviation  of  the 
temperature  near  the  laminar  sublayer,  as  found  using  heat  transfer  theory. 
This  is  discussed  in  Appendix  C.  We  see  from  Eq.  (33),  again  near  the  sub- 
layer, that  the  temperature  standard  deviation  is 

/ 

= 0.  027(T^  - TJ  M,  (35) 

with  a nominal  value  used  for  the  turbulence  intensity  and  T„/(Z')AZ  replaced 
by  T^  - T^.  This  result  is  approximately  a few  percent  of  the  temperature 
contrast.  The  Mach  number  dependence  is  very  weak,  probably  quite  diffi- 
cult to  find  experimentally.  The  dependence  on  distance  downstream  from 
the  leading  edge  likewise  is  very  weak. 

The  normalized  temperature  correlation  is  given  in  Eq.  (34).  It  seems 
plausible  thait  the  x -dependence  of  the  correlation,  not  found  from  the  present 

analysis,  will  be  essentially  similar  to  the  dependence  in  the  y-direction. 

2 2 2 

We  propose  replacing  T)  + f in  Eq.  (34)  by  a defined  by 

<j  = (x  +y  +z  )/M  (36) 

where,  recall,  M is  the  local  scale  of  the  turbulence,  approximately  the 
distance  from  the  wall  when  close  to  the  wall,  and  6 is  the  displacement 
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thickness  of  the  boundary  layer  as  one  moves  out  into  that  layer.  Then  the 
normalized  correlation  may  be  written  following  Eq.  (34) 


T(r\  r'*) 

[tV)]^ 


= Ij(a^  + C%(a^) 


(37) 


To  calculate  these  integrals  we  need  to  know  the  normalized  longitudinal  velo- 
city correlation  f.  This  function,  for  the  large  Reynolds  number  flows  of 
interest  here,  could  be  taken  from  experiments;  however,  we  must  differen- 
tiate the  function  once  in  constructing  the  required  integrals.  That  makes  the 
process  somewhat  more  difficult.  Therefore,  we  have  chosen  to  take  the 
function  f from  recent  large  Reynolds  number  calculations.  The  longitudinal 
correlation  function  of  Hogge  and  Meecham  has  proper  continuity  and  dif- 
ferentiation is  not  a problem.  The  results  of  the  computation  of  the  two  in- 
tegrals and  are  presented  in  Fig.  Z.  Specifically,  for  the  two  measuring 
points  at  the  same  height  z'  from  the  plane,  the  correlation  function  is  just 
I^.  It  is  interesting  to  note  that  there  is  a small  negative  lobe  in  normalized 
value  of  about  -0.2  at  a distance  in  the  x-y  plane  of  Z to  3 local  correlation 
lengths.  This  arises  from  the  known  negative  lobe  of  the  transverse  velocity 

correlation,  on  which  this  integral  is  dependent.  If  we  look  at  the  correla- 

2 2 2 

tion  in  the  vertical  direction,  that  is  x = y = 0,  we  need  Ij(C  ) + C ■^S*^*^* 

it  appears  that  a negative  correlation  will  develop,  according  to  our  results, 
though  it  is  not  large. 

The  density  fluctuation,  which  is  needed  for  optical  propagation  through 
the  heated  boundary  layer,  may  be  obtained  from  the  ideeil  gas  law.  Recall 


I 

I 
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that  the  variations  are  at  constant  pressure,  therefore,  it  follows  that 

p/p(z')  = -T/T(z')-  (38) 

Here,  p(z^)  is  the  local  average  value  of  the  density  and  p is  the  density  fluc- 
tuation. It  is  easily  seen  that  the  normalized  density  correlation  function  is 
given  by 


Again,  we  assume  that  the  intensity  of  the  temperature  fluctuation  ^nd  mean 
values  of  the  temperature  vary  but  slightly  in  going  from  t!  to  7,"  in  regions 
of  interest.  The  normalized  density  correlation  is  the  same  as  the  normalized 
temperature  correlation.  Further,  the  rms  value  of  the  pressure  and  tempera- 
ture are  related  from  Eq.  (38)  by 


p'(z')/m*')  = T"(z')/T(z') 


(40) 
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APPENDIX  A.  THE  EVALUATION  OF  THE  LAGRANGIAN 
INTEGRAL,  Eq.  (16) 

It  is  seen  in  th^  body  of  the  text  that  we  attempt  to  evaluate  a Lagrangian 
integral,  that  is,  an  integral  of  the  velocity  fluctuation  following  the  flow.  The 
exact  form  of  that  integral  can  be  written 


T(r)  = >T^(z) 


where 


/•'^D  r /•'^D  ^ /•'^D 

J w jx  (u  + u)df,  y - y vdt, 

/•^D  ^ 1 ^ 

z - / wdt;  t - Y I dt, 

- -f  J 

u = u z - # wdtj 


r r^-D  ^ /•^D  1 

u = u|x  - y (XT  + u)  d?,  y - J vdt,  ^ - J wdE’;t  - ?J  (A2) 


and  similarly  for  v and  w.  Such  a Lagrangian  integral  is  clearly  extremely 
complicated.  There  is  a continued  implicit  dependence  upon  the  Eulerian 
space  variables  x,  y,  z,  and  the  time  t.  In  Eq.  (16),  we  broke  off  that  ex- 
tended implicit  dependence  by  neglecting,  in  the  arguments,  the  velocity 
fluctuation  u. 

We  attempt  here  to  find  the  qualitative  effect  of  that  approximation.  The 
most  important  such  effect  can  be  seen  in  Eq.  (A2).  By  neglecting  w in  the 
argument  of  the  function  XT,  we  neglect  changes  in  the  vertical  position  which 
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may  cause  large  changes  in  the  value  of  the  average  flow  velocity  experienced. 

Referring  to  Eq.  (Al)  it  is  seen  that  these  in  turn  cause  relatively  large  changes 
in  the  x argument  of  the  function  w.  This  introduces  a random  fluctuation  in 

■] 

the  x-direction  which  has  been  neglected.  A consequence  of  that  neglect  is 
that  the  correlation  of  the  temperature  fluctuation  in  the  x-direction  assumes, 
incorrectly,  large  scales  in  that  direction  as  was  seen  in  the  body.  Further, 
as  a result  of  our  inability  to  deal  exactly  with  the  complicated  integration 
represented  by  Eq.  (Al)  to  (A2),  there  will  be  errors  in  the  amplitude  of  the 
predicted  fluctuation,  mostly  on  the  large  side,  since  we  have  suppressed 
some  of  the  fluctuating  character  of  the  integral  in  (Al).  As  a consequence, 
it  was  seen  in  the  text  that  it  was  necessary  to  introduce  an  empiriced  coef- 
ficient in  order  to  match  the  known  temperature  fluctuations  near  the  laminar 
sublayer. 


i 

» 


-30- 


i 


There  is  a question  concerning  whether  the  form  Eq.  (22)  is  a proper 


one  for  an  inhomogeneous,  incompressible  velocity  field.  We  know  that,  for 
homogeneous  turbulence,  Eq.  (23)  will  yield  such  an  incompressible  field,  but 
the  weighting  fvinctions  appearing  in  Eq.  (22)  complicate  the  discussion.  Con- 
sider the  relation 


R,.(r',  r')  = sl;;-  *ij  ’ ^ ^jw'(z' )w'(z'' )M^z', 

X /'%'V)dp' 

•4-k 


tT) 


(Bl) 


where  we  use  the  sunrunation  convention:  repeated  indices  should  be  summed 

from  one  to  three.  For  incompressible  flow  the  divergence  of  the  correlation 

3 

function  must  vanish;  thus 


^ Ri.(r',  r")  = ^ Rij(r'.  rO  = 0 


(B2) 


It  is  readily  verified  that  Eq.  (Bl)  has  this  property.  Furthermore,  it  is 
easily  seen  that,  for  statistically  homogeneous  turbulence,  where  y/  and  M 
are  constant,  Eq.  (Bl)  reduces  to  the  homogeneous  form  given  in  Eq.  (23). 
Finally  for  the  z-component  correlation  we  Bnd  the  form  proposed  in  Eq.  (23). 
Thus,  the  correlation  Eq.  (Bl)  is  incompressible,  reduces  to  the  homogeneous 
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correlation,  and  does  yield  the  simple  form  proposed  in  Eq.  (23).  It  is  seen 
that  correlations  other  than  have  coe^icients  involving  the  turbulence - 
intensity  weighting  function,  V,  in  a more  complicated  way. 


APPENDIX  C.  CALCULATION  OF  THE  CONVECTED  HEAT 


FLUX  q AND  T"(fl^) 


The  beaic  equation  for  the  calculation  of  the  convected  heat  flux  is 


q(«,)=»(T.„-TJ  . torT.„>T„ 

= NT,  - T.)  . for  T.,  < T, 


essentially  Newton's  law  of  cooling.  Here  T^^  is  the  temperature  which 


would  be  attained  by  a thermally  insulated,  adiabatic,  wall. 


result  of  the 


heat  generation  within  the  laminar  sublayer.  The  constant  of  proportionality 
h is  obtained  from  the  empirical  relation 


h = Pr*'^^^p*u^C*  0.0288(Re*)"*^® 


A discus sioei  of  this  type  of  heat  transfer  calculation  may  be  found  in  Ref.  4. 
The  usual  calculation  is  for  the  heat  transfer  at  the  wall  (what  is  needed  in 
previous  applications).  We  use  here  the  same  coefficient  h,  since  the  convec- 
tive process  needed  is  in  the  vicinity  of  4^.  This  transfer  coefficient  applies 
to  turbulent  boundary  layers t 5 X 10^,  the  case  here.  It  will  be  seen 

that  the  ratio  h/C^p  (which  we  need)  is  but  weakly  dependent  on  tite  parameters 
of  the  calculation. 

. t 

The  starred  quantities  in  Eq.  (CZ)  are  all  referred  to  a reference  tem- 


perature giSen  by 


T*  ■ T4  + 0.§(T^  - 


TJ  + 0.Z2(T^^  - T.) 
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The  Prandtl  number  Pr  is,  for  air,  0.71,  and  the  Reynolds  ntonber  is  defined 


Re*  = 


o> 

V* 


(C4) 


The  stagnation  temperature  T^  is  defined  by 


(C5) 


where  y is  the  ratio  of  specific  heats.  A recovery  factor  giving  the  effect 
upon  the  heat  transfer  of  the  process  of  bringing  the  high  speed  fluid  to  rest 
(through  viscovis  action)  at  the  wall  is 


r 


= Pr 


1/3 


(C6) 


The  procedure  for  tiie  calculation  of  q is:  first,  find  the  recovery  fac> 
tor  r and  the  stagnation  temperature  T^;  then  determine  the  reference 

tennperature  T*,  the  Newton's  law  constant  h,  and  finally  use  Eq.  (Cl)  to 
find  q. 

It  should  be  noted  that,  for  flows  such  as  this,  tiie  pressure  remains 
essentially  constant  throughout  the  boundary  layer,  for  fiiis  calculation. 

For  our  purpose,  we  can  simplify;  divide  by  C p(T  ) and  separate  out 

p •W 

the  M.b  u^/e.  dependence. 


Pr'^^^M^®  0.0288(a,/v*)'‘^®a,x'^^® 


(C7) 


-34- 


where  x is  measured  from  the  point  of  initiation  of  the  boundary  layer  and  we 
use  the  fact  that  Pr  is  nearly  const2Uit. 


Consider  an  example:  a jet  plane  flies  at  an  altitude  of  28  Kft,  at 

M = 0.9  with  T = 72*F.  At  this  altitude  T = -40*F,  p =0.30  atmospheres; 

we  want  to  know  h at  x = 20  ft  (we  use,  as  usual,  English  units).  From  the 

relations  given  above  we  find  the  intermediate  values  a^  = 1006  ft/sec, 

T = 488*R,  r = 0.892,  T = 481*F,  T*  = 489*R,  Re*  = 4.65  X 10^,  v*  = 
o aw  X 

4.33  X 10'^  ft^/sec.  We  have,  for  Eq.  (C7),  converting  to  MKS 


h/C  p = 1.52X  lO’^M^^a 

pr  o»  » 


= 1.40  X lO'^a.x*^''^,  in  MKS 


We  see  T > T , so  from  Eq.  (Cl) 
w aw  n \ f 

q/CpP  = 1.52X  10-^My®a„x'^^^T^  - T„)*C  - m/s 
= 18.  6*C  - m/s 


(C9) 


Consider  a second  example,  a nominal  wind  tunnel  experiment,  take 

= 0*F,  = 0.8,  p^  = atmospheric  pressure,  = 70*F.  We  find 

T = 519’R,  T = 513*R,  a = 1053  ft/sec,  T*  = 507’R,  v*  = 1.50X 
o aw  ® 

♦ 7 

sec;  choose  x = 3.  7 ft. , then.  Re  = 2.08  X 10  . We  have 

X 

h/C  p = 1.22  X lO’^M^^a  MKS  (CIO) 

n ' QD  OD 
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We  see  T > T ; thus, 
w aw 

q/CpP  = 12.  4®C  - m/s  (Cll) 

4/5  -1/5 

We  note  from  Eqs.  (C8)  and  (CIO)  that  the  coefficient  of  a^x 

chcUiges  little  in  these  quite  different  applications;  we  might  for  many  appli- 
cations of  the  above  type  take  it  as  constant.  Thus,  for  such  application  one 
can  use 

q/CpP  = 1.3X  - T„)“C  - m/s  (C12) 

in  MKS  (and  *C). 

Note  also  that 


= <wT>~  w'T' 


(C13) 


« 


if  <wT(£  = 0)>  ~ w^T^.  There  is  the  possibility  that  the  average  is  much  smaller 
than  that  product;  in  statistically  homogeneous  and  isotropic,  incompressible 
turbulence  <wT(r^)>  = 0,  the  whole  correlation  vanishes  though,  of  course, 
v/  and  T'  do  not.  All  things  considered,  however,  Eq.  (Cl 3)  seems  plausible. 
So  using  w^  /u^  0.  05,  we  obtain 


'!'(«) 


1.3  X 10-V/V(x/x„)'*^® 


■0'8u_ 


(T  - T„) 
' ^ 00 ' 


(C14) 


= 0.  026  M.(x/x^)'*''^T^  - T.)  , x^  = 1 m 


I 
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The  maximum  temperature  flucttiation  is  a few  percent  of  the  temperature 
contrast,  as  one  might  have  assumed.  These  results  are  clearly  not  valid 
for  small  Mach  numbers  or  small  x. 


THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  t.pccies  in  rocket  piumes.  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  .hemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory;  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution:  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue- induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields:  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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